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Abstract. We study the collective two-channel dissociation dynamics of a molecular 
Bose-Einstein condensate into bosonic fragments under tight harmonic confinement. 
Bose-stimulated dissociation in either channel can only take place provided that the 
respective trap size U for the fragments is large with respect to the healing length Q 
of the atom-molecule resonance. Thus, even when both channels are equally coupled, 
differences in mass or in dynamical polarizability enable the control of the reaction 
outcome by variation of the trap frequency. In particular, if l\ > Q\ and I2 < C2, 
only the first channel will be amplified. This behavior demonstrate a unique feature of 
'supcrchemistry' wherein a chemical reaction may be controlled by the manipulation 
of the reaction vessel. 
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1. Introduction 

Quantum degenerate bose gases of dimers consisting both fermionic [U [21 El El E] and 
bosonic 0, El [8] constituent atoms, have recently been realized in Feshbach resonance 
P Uni Ell H2] and photoassociation [131 HH UHl HE1 HEl UHl EQl EH [22] setups. First 
evidence for boson pairing was obtained in seminal experiments in a BEC of 85 Rb atoms 
[23] showing Ramsey-like oscillations indicative of the coherent coexistence of molecular 
and atomic condensates. This experiment was soon followed by the production of much 
larger and well-defined molecular condensates of 133 Cs2 [6], 87 Rb2 [7], and 23 Na2 [8], as 
well as ultracold molecular gases and molecular BECs with fermionic constituent atoms, 
such as 40 K 2 [1] and 6 Li 2 [21 El El E] . BECs of molecules consisting fermionic atoms,turn 
out to be more stable due to the Pauli blocking of collisional relaxation [24j . 

These novel entities open the way for the corroboration of theoretical proposals 
for, so called, superchemistry [23, EEJ [23, EHl EHl [30]. The underlying dynamics of 
atom-molecule conversion in quantum gases, is collective in essence and is therefore 
strongly affected by the quantum statistics of atoms and molecules. For the purely 



Confinement controlled dissociation of a molecular Bose-Einstein condensate 2 

bosonic system, bose enhancement takes place for atoms and molecules, to result in a 
nearly exact analogy [251 ED E21 E31 El] with parametric processes in optics, including 
coherent matter wave amplification, phase conjugation, and quantum state squeezing. 
Association of atoms to form a condensate of dimers is analogous to second-harmonic 
generation and the dissociation of a molecular condensate corresponds to parametric 
downconversion. The latter pair-production process is known to exhibit a modulational 
instability, resulting in from the exponential amplification of spontaneously emitted 
photon-pairs [351 ESI E3 EE] • 

The existence of similar parametric gain in molecular dissociation was noted early 
on in the theoretical study of atom- molecule condensates [25J E2 ETJ EHl EH 126] . 
Moreover, since driving the system in the vicinity of a dynamical instability generates 
gain in quantum noise [39], there are significant deviations from the classical- field 
predictions of the Gross-Pitaevskii (GP) mean-field theory [321 Hd HD H21 03]. The 
rapid growth of correlations can be used for the generation of pair-correlated or number- 
squeezed atomic beams [HI 05] using either spontaneously emitted pairs [B] or a 
classical atomic (signal) seed [15] , in analogy to the use of parametric downconversion 
to generate squeezed light in quantum optics 

Recently, we have found that stimulated molecular dissociation can be arrested, 
provided that atoms are confined to move within an harmonic trap of frequency 
ujq, whose size Iq — a/ hf (mono) is smaller than the resonance healing length ( = 
h(mg)~ l l 2 n~ l l 4: [17], where m,g, and n are the atomic mass, the atom-molecule coupling 
strength, and the average molecular density, respectively. Thus, the large coherence 
length of the molecular BEC on the one hand, and the tight confinement available 
today by magnetic and optical traps on the other hand, conjure to produce collective 
chemistry which not only depends on particle statistics, but also on the size and 
shape of the 'reaction vessel' itself. Here we show that for two-channel dissociation 
of a molecular condensate of triatomic ABC molecules, the critical dependence on the 
trapping potential can be used to control the dissociation outcome by mere manipulation 
of the trap parameters, even when the two rearrangement channels are equally coupled 
to the molecular BEC. Thus for sufficiently tight traps, neither channel is amplified. 
At the other extreme, for trap size larger than the resonance healing-length of both 
channels, both channels will be amplified. In between these regimes, lies a window 
where only one channel attains exponential gain, resulting in nearly perfect selectivity. 

The exponential gain features of molecular BEC dissociation and the possibility to 
arrest the dissociation of a homonuclear diatomic quantum gas by tight confinement 
are presented in Section [2l The extension to the case of heteronuclear molecules where 
dissociation fragments have different masses and optical polarizabilities, resulting in 
differences in trap-size and healing length, is provided in Section El In Section 0] we 
present the analysis for the two-channel dissociation of a triatomic quantum gas into 
bosonic constituents and corroborate the analytical prediction by numerical examples. 
The values of experimental parameters required for the observation of the predicted 
phenomena, are briefly discussed in Section [51 while conclusions and future directions 
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are presented in Section 



2. Exponential gain and stabilization by confinement 



In this section we provide a brief review of the results of Ref. [U] which serve as the 
foundation for confinement control. We consider the dissociation of a molecular BEC 
constituting dimers of identical bosonic atoms. The atom-molecule Hamiltonian in the 
rotating wave approximation is, 
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where ^(r, t) and 0(r, t) are atomic and molecular bose field operators, m is the atomic 
mass, A is the detuning from atom-molecule resonance (i.e. the interaction-free energy 
difference between the threshold of the open channel and the molecular bound state 
in the closed channel of a Feshbach-coupled system, or the detuning from two-photon 
resonance in a photodissociation scheme) and V^(r) = muoQ 2 r 2 /2 is an isotropic optical 
trap potential with frequency u Q . Atoms and molecules see essentially the same trap 
frequency since molecules have twice the mass and twice the optical polarizability [48J. 
The molecule-molecule interaction strength is U m and g is the atom-molecule Feshbach 
or optical coupling frequency, assumed without loss of generality to be real. Interactions 
between atoms are neglected since we are interested in the initial dissociation dynamics 
when atomic densities are small. 

The exponential gain features of molecular BEC dissociation are readily captured 
within an undepleted pump approximation. This amounts to linearization of the atom- 
molecule Hamiltonian (OQ) about a constant molecular BEC wavefunction 0(r, t) — > 
(f)(r)e~ ifJ,t / h , obeying the Gross-Pitaevskii time-independent equation 

2 - J .+ [2V(r) + t/ m |0| 2 ]0 = /i0. (2) 



4m 

The dynamical equations for atomic field operators \l/ 
read, 
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where 5 = A — /i/2 is an effective detuning including molecular zero point kinetic energy 
and molecular mean- field shift. 

This linearization procedure (which is used extensively in three- and four-wave 
mixing in optics [351 1361 [371 138] ) is inherently limited to short-time dynamics where 
molecular condensate depletion as well as fluctuations of the molecular field are small. 
Its validity regime is set by the standard parametric-approximation limits: gt — > 0, 
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— * oo while gcj)t = constant [3B] • Within these limits, it constitutes a reliable stability 
analysis towards stimulated dissociation. Molecular condensate stability is determined 
by the characteristic eigenvalues A of the set (J3j)- (SI) - For a uniform gas with V(r) = 0, 
= \fn, and \i = U m n with n being the molecular condensate density, the system 
is readily diagonalized in momentum space ?p(r,t) = a>q{t) exp(iq ■ r), to give the 
dispersion relation, 



xj \e q + 5) 2 - n\g\ 



(5) 



where e q = [hq) 2 /2m is the free-particle dispersion. Thus for a uniform gas with 
5 < y/n\g\, there are always unstable resonant modes for which e q + 5 < gy/n, resulting 
in complex characteristic frequencies and exponential gain. The dynamics of amplified 
modes is described by well-known solutions, given in terms of the hyperbolic functions 
cosh(A g t) and sinh(A (? t), which can be found in Refs. [321 HH HH1 HE] ■ F° r zero effective 
detuning 5 = 0, the unstable modes are simply low-energy, long- wavelength excitations. 

The stability analysis for a nonuniform molecular gas can be gleaned from the 
Bogoliubov stability analysis of an attractively interacting atomic BEC with the 
Hamiltonian, 



H 



—V 1 + V(v) + |^ 



ipdr . 



(6) 



The Bogoliubov equations are obtained from the usual separation ip(r,t) = ip(r) + 
5ip(r,t), where ip(r) is the c-number solution of the time-independent GP equation 
[-|^V 2 + V(r) + U\^\ 2 ]i; = flip and 5ip(r,t) is an operator correction. Linearization, 

obtained by retaining up to quadratic terms in Sip, Sip and rotating 5^1 = S^e 1 ^^, 



5V = Sip e 



-i/it/h 



gives the familiar forms, 
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Equations are similar to Eqs. (j3J- ([!]), except for the replacement of the classical 

molecular pump by the atom density ip 2 and the addition of the diagonal mean-field 
term 2U\ip\ 2 corresponding to normal pairing of the $ ^ tpip interaction, which does not 
exist for the Feshbach coupling (j^ipip + H.c. The Bogoliubov transformation S^f(r, t) = 

£ q ( 



a q w q (rje 



6t* v * 



>)e ?Aqt ) seeks to diagonalize the set ([7j)-(jH]) similarly to the 
diagonalization of Eqs. ([3])-(Jl]). For a uniform bose gas in 3D, the additional term results 
in the modification of the dispersion relation to the famous Bogoliubov dispersion 



frA q = ^J(e q + nU) 2 -(nU) 2 = yj e q (e q + 2nU). 



(9) 



For an attractively interacting gas with U = 4nh 2 a s /m < (a s being the s- 
wave scattering length of the atoms), long wavelength excitations with e q < \2nU\ 
corresponding to phonon-like excitations for equal-magnitude but positive U, attain 
complex frequencies and are thus unstable, leading to the collapse of the condensate 
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[1HJ EH EI]- However, for a BEC subject to a confining potential, collapse can be 
arrested provided that the size of the trap Iq is short with respect to the healing length 
f , i.e. 



h < / h2 ~£ (10) 
mcjo y 2mn\U\ 

where n denotes the average condensate density. Thus, provided that the spacing 
between trap levels TvjJq is larger than the average interaction energy |2nC/|, the unstable 
modes are precluded from the trap and stable condensates of attractively interacting 
atoms exist up to a critical density for a given trap- frequency [52} E3] . The maximum 
number of particles allowed for a given trap, is obtained from writing the trap frequency 
and the interaction energy in terms of the trap size Iq, 

H 2 N N 
f)ujQ = ? , Un = U—^U ~, 11 

™Iq 2 V ( 4 /3)tt (V2i ) 

where V denotes the trap volume, corresponding to the volume of a 3D spherical box 
whose radius is equal to the characteristic size of the interaction- free ground state V21q. 
This gives the familiar stability condition 

N<Nc ^^L. fJtL, (12) 
3 2m\U\ 3 \a s \ 

stating that the critical number of particles in the attractively interacting BEC 
corresponds to the number of s-wave scattering lengths contained within the 
characteristic trap size. For 7 Li in harmonic traps of order 100Hz, l is of the order 
of a few fim, compared with the ~ 30 bohr scattering length, giving critical numbers 
Nq of several hundreds to few thousands particles [53]. Interestingly, experimental 
measurements give Nc\a s \/lo = 0.46 ±0.06 [54]- This value is in an excellent agreement 
with the prediction above of y/2/3 = 0.47, despite the fact that the theoretical value is 
based on a crude approximation for the density. However, due to the rather arbitrary 
choice of trap radius (taken here to be the point at which the ground-state order 
parameter drops to 1/e of its maximal volume), the agreement appears to be incidental 
rather than profound. 

Returning to molecular dissociation, a similar threshold for instability is obtained, 
despite the differences between the gapless Bogoliubov dispersion relation and the 
'imaginary-gapped' molecular dissociation spectrum (jHJ). Provided that the trap size is 
smaller than the 'resonance healing length' (, i.e. 



k = J—<J-^m = C, (13) 
V mujQ y mgn 1 '^ 

or alternatively if frwQ > g\ffi, there will be no unstable modes towards dissociation, 
even when 5 = (i.e. on-resonance coupling). The resonance healing length is the 
characteristic lengthscale of Eqs. ([3]) and (j4j), resulting in from the interplay of kinetic 
and resonance energies. Thus, when this lengthscale exceeds the system size imposed 
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Figure 1. Confinement effects on homonuclcar dissociation: (a) Real and imaginary 
parts of the characteristic dissociation frequencies as a function of atom-molecule 
coupling strength, for an harmonic trap with frequency ojq. The coupling is resonant 
5 = 0, and the molecular condensate is nonintcracting, U m = 0. Complex frequencies 
appear for g\ffi > ha>o (b) Total atomic fraction obtained in a mean-field dynamical 
calculation, triggered by Gaussian zero-average noise, as a function of rescaled time 
for g^/E/hojo = 0.9 (dashed red line) and for g^/E/hujQ = 1.0 (solid blue line), with 
S = U m = (c) Atomic density profiles from the same dynamical simulation at t = 30. 
Amplification is clearly attained above the expected threshold. 



by tight confinement, the difference in zero-point motion between bound molecules 
and pairs of dissociating atoms, inhibits the dissociation process. This effect is a 
characteristic of 'Superchemistry' [26] in that the outcome of a coherent chemical process 
depends due to its collective nature, not only on the statistics of the constituent atoms, 
but also on the size and geometry of the 'container' in which it is carried out. It results 
in from the fact that the entire gas behaves essentially as a macroscopic single-particle 
with a relatively large coherence length, and from the possibility to introduce confining 
potentials on the same /xm lengthscale. 

In Fig. UK we plot the characteristic frequencies of Eqs. ©-(H]) in ID with 
^ = U m = and a finite-size undepleted pump obtained from Eq. (j2J), by expanding 
0(r) = Y^j { PjXj( r )i ^( r :t) — Ylij a j(.t)Xj( T ) i n the basis of harmonic-oscillator solutions 
Xj( r ): an d diagonalizing the resulting linear set of equations for the dynamics of Oj(i) 
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and dj(t), 




with a = (S , Si, ...), a 1 " = (S , Sj, fit, ...), = u 5 i:j (j + 1/2), = 
(,9/h)2~2i i Pi{XjXi\Xi)- Identical eigenfrequency plots were obtained from a discrete 
Fourier transform diagonalization. Looking at the eigenfrequencies dependence on the 
atom-molecule conversion frequency g\/~fijh, it is evident that complex frequencies, 
corresponding to an initially exponential molecular gain, are only obtained if the 
coupling frequency is larger than the trap frequency. Thus, for g\ffijfvjjQ < (3 where f3 is 
a geometrical factor of order one, the molecular condensate is stabilized by atomic zero 
point motion. The ensuing dynamics, including molecular BEC depletion, is depicted 
in Fig. [Tb where atomic fraction is plotted as a function of time for g^/E just below 
the predicted amplification threshold (dashed red line) of g c ~ 0.91huj/ \ffi and slightly 
above it (solid blue line). The equations were solved in a mean- field approach, simulating 
initial spontaneous dissociation by a complex Gaussian white noise with zero mean, 
corresponding to atomic quantum fluctuations in the truncated Wigner representation 
[551 156] . This technique is not as complete as full stochastic calculations [SI US], but 
it is sufficient to accurately predict the onset of instability. For g < g c all frequencies 
are real and there is no observed dynamical gain, whereas for g just above the critical 
value, molecular amplification takes place. The corresponding atomic density profiles 
at time u t = 30 are plotted in Fig. [Tfc. 



3. Dissociation of a heteronuclear molecular BEC 



Having established that molecular BEC dissociation depends critically on trap 
parameters and that the stimulated process can be arrested altogether by sufficiently 
tight confinement, we turn to the case of heteronuclear diatomic molecules. The main 
difference from the homonuclear case is that due to differences in mass and dynamical 
polarizability, the dissociation fragments see effectively different trapping potentials and 
have different healing lengths. Thus we need to reformulate the amplification conditions 
accordingly. 

We consider the dissociation reaction AB — > A + B, where diatomic AB molecules 
are fragmented into their A and B constituent atoms, in an off-resonant optical trap. 
We denote the atomic masses by m^, mg and the atomic dynamical polarizabilities by 
a a and respectively. The atomic trap frequencies experienced by the atoms are 
proportional to the square root of the polarizability to mass ratio, 

^A,B °c • (15) 

m A ,B 

Since the polarizability of the loosely bound AB molecules is approximately a a + «b, 
and their mass is niA + m B , the molecular trapping potential will be given to a good 
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Figure 2. Elements of the coupling matrix <£>„„/ with U m fi/hujAB = 10, for itia = mB 
(a) and for tub = 2m a (b). The even parity of the molecular condensate wavefunction 
forbids coupling between trap states with opposite parity. Diagonal elements dominate 
the coupling matrix. 
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reducing as required, to the atomic trap frequency when oj a = uj b = 

Given these trap frequencies, the dynamical equations for the atomic field operators 
^ a,Bi linearized about the molecular order parameter 0(r)e~ JM */ ft obtained from the 
stationary GP equation, 
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In what follows we shall assume 5 = 0. Rescaling length as r ^ t/Iab with 
Iab = \fh/[(jriA + ^b)^ab\, normalizing $ = 4>/\fN (so that f |$(r)| 2 <ir = 1), and 
defining the dimensionless time r = u>ABt, Eqs. (|17p - (|19|) transform to 
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where M = + is the molecular mass and n = N/lAB d is the average molecular 
density with d denoting the number of dimensions. The characteristic frequencies of the 
set (l2Tj) - (l22]) are obtained by expanding the atomic field operators in terms of eigenmodes 
of the respective atomic traps, e.g. in ID, 
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reduces the dynamical equations into a coupled set of equations for the mode 
annihilation operators a n and b n , 
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where <& nn / = J dxx^(^)$(^)Xn'( x )- The coupling matrix elements $ nn / are plotted in 
Fig. [2] for 1 < n,n' < 10. Since the molecular condensate wavefunction possesses even 
parity, only trap modes with equal parity are coupled by the Feshbach interaction term. 
Moreover, because the molecular density is close to the ground state of the atomic trap, 
the diagonal elements dominate the coupling matrix. Neglecting all off-diagonal terms 
(this resembles the coupling of the the modes q and — q in the uniform case), we arrive 
at the characteristic frequencies, 

(27) 



±n 



2 L 



n 



0,1,2, 



where e^' B = fruJA,B( n + 1/2) are the respective harmonic oscillator energy levels. Thus, 
if the effect of off diagonal $ n / n elements can be neglected, amplification will be attained 
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Figure 3. Confinement effects on heteronuclear dissociation: (a) Real and imaginary 
parts of the first four characteristic dissociation frequencies, with U m fi/hu>AB = 10, for 
vtia/cia = wib/olb (solid blue line) and for ttiaIola = TOb/2q;b. The critical value for 
stability to dissociation depends on the ratio ttibola /niA&B- (b) Dependence of critical 
coupling strength on mass ratio m^/mA for U m n/hiUAB = 10 and a a = olb- The exact 
value obtained by direct diagonalization (solid green line, circles) is compared with the 
approximate critical coupling from Eq. (|28| (dashed black line). 



only if \fng(^ nn > H{uja + Wb)(ti + 1/2) for some n. Since $ nn < 1 for all n, we obtain 
the amplification condition, 

IV^I > h ±A±^A . (28) 

In Fig. [3^ we plot the first four characteristic frequency obtained by full 
diagonalization including all off diagonal <3> nri / terms, as a function of the coupling 
strength gy/n, for niA — vtib (solid blue line) and for = 2m a (red dashed line). 
The optical polarizabilities were assumed for simplicity, to be identical a a = olb- 
The resulting eigenvalues turn out to be almost identical to the approximate formula 
(127)) . For equal masses, we have uja = ojb = ujab, so that amplification is attained 
for gy/E > Jtljab/^oo = Zfo^AB- We note that the factor two difference from the 
homonuclear case results in from the distinguishability of the dissociation fragments. 
A similar factor of two separates the nondegenerate parametric amplifier from its 
degenerate counterpart. For m# = 2m a and a a = cxb we have uja = V^oob = \/3/2ujab, 
so that uja + ujb = ^Kj±^) Uab ~ 2.1ujab- Consequently, complex frequencies begin 
to appear for g\/fi ~ 2.1Hujab, as shown in Fig. [3k In Fig [3b the critical coupling 
strength for stimulated dissociation, is plotted as a function of the mass ratio ms/m^. 
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Figure 4. Real and imaginary parts of the first eight characteristic dissociation 
frequencies for two-channel dissociation with thia/oia = tub/olb = mc/'Z&c an d 
U m n/hujAB = 10, as a function of the coupling parameter g\ffijhiOABC- Solid blue 
lines denote the frequencies associated with A+BC channel whereas dashed green lines 
correspond to the AB + C channel. Vertical dotted lines mark the critical coupling 
strength of each channel. In between these lines, amplification will be attained only in 
the AB + C route. 



The solid green line is obtained by direct diagonalization, whereas the dashed black 
line corresponds to the approximate amplification condition (128]) . evidently giving an 
excellent estimate of the amplification threshold. Numerical simulations of dynamics, 
carried out just below and just above the anticipated critical coupling, clearly show the 
transition from small oscillations and an initial exponential gain. Thus the outcome of 
dissociation will depend on the mass ratio between the dissociation fragments. 



4. Control of two-channel dissociation by confinement 

The possibility to control multichannel dissociation through confinement is a 
straightforward extension of the previous sections. Consider a quantum gas of triatomic 
ABC molecules which may undergo dissociation into either A + BC or into AB + C [28J. 
This is a second- quantized version of the familiar coherent control paradigm. Given that 
the constituent atoms are subject to trap frequencies lja,b,c which are proportional to 
^/ a A,B,c/ m A,B,c respectively, the trap frequencies for the BC and AB fragments can 
be approximately given as u 2 BC = mB ^+^T c and u 2 AB = mA ^_ZT B ' res P ectivel y- 
The trap frequency for ABC molecules is uj ARr = mAU} A +m BU B +m c u c ^ jyj agg differences 

r ~l j ah^j mA -\-m,B +mc 
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between the constituent atoms translate to differences between the m^/fmg +mc) and 
the mc/ (t^a + ^b) fragment mass ratios, resulting in different amplification thresholds. 
Consequently it is possible to control the dissociation by variation of the trapping field. 

Rescaling as before r — » r/l ABC with Iabc — \fhj (Mluabc), t — > r = u>ABct, 
(p — > $ = (f>/y/N, the two-channel linearized equations (assuming Ai = A2 = yu/2 with 
Ai ; 2 being the respective A + BC and AB + C channel detuning), are a simple extension 
ofEqs. (EDD-d22]), 

; $ = (29) 
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Here gl and #2 denote the atom molecule coupling strength in the A + and AB + C 
channel, respectively. Characteristic frequencies are again found by expanding all the field 
operators of the fragments, in terms of the respective trap modes and diagonalizing the 
resulting set of equations for the mode annihilation and creation operators. Since the 
two channels are only coupled through the mutual pump field, the eigenvalue problems 
separate and the approximate eigenfrequencies are the same as in Eq. ([27)1 . 
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where $ 1>rm # = / dxx A (x)^(x)x B F ! (x), $2, n «' = / dxx AB {x)${x)Xn>{ x ) and Xn( x ) are 
the respective oscillator wavefunctions for species a. Consequently, the amplification 
condition for each channel is roughly, 



00 



> 



h(uj A + ujbc] 



n\g 2 ® 2 , 



on 
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(35) 



2 ' ' 1 2 

By judicious choice of the confining potential it is thus possible to find regimes where 
gain is attained only in one of the two dissociation channels, even when g\ = g 2 and both 
channels are equally coupled to the molecular BEC. As an example, we shall consider an 
arbitrary case where m^/a^ = ms/ctB = mc/2c*c, Q\ = gi = g, an d U m n/hujABC = 10, 
corresponding to [i — 3A1Hujabc- In Fig- BJ the exact characteristic frequencies h\±„ BC 
(solid blue line) and \± B +c (dashed green line) are plotted as a function of the coupling 
strength. The critical values for amplification are g^, BC \f^,jhu>ABC — 2.1960 and 
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Figure 5. Fraction of dissociated atoms in the A+BC channel (solid blue line) and the 
AB + C channel (dashed green line) as a function of rescaled time r = u>ABct- Atomic 
masses, optical polarizabilitics, and the molecular interaction strength are the same as 
in Fig. [4j Coupling strengths arc: (a) g\ffi / 'fuv abc — 1-9, (b) gy/E/fkuABC =2.1, (c) 
g\f^/ftwABC = 2.3. In (d)-(f) we retain the coupling strength, length and time units 
of (b) and multiply all trap frequencies by a 'widening factor' w of 0.8 (d), 1.0 (e), and 
1.2 (f). The channel population distribution depends on trap size. 



9ab+c^ l /^ JjJ abc — 1-9960, in very good agreement with the approximate conditions 
(I35p . In between these critical coupling values, there lies a window where exponential 
gain will be attained only in the AB + C channel. This prediction is supported by 
the numerical simulations of dynamics (including pump depletion) for various values 
of gy/E/hu) abc, shown in Figs. Eh through [5b. When g\ffijhjjJABC = 1-9 (Fig. [5k) 
the coupling is subcritical and both channels are stabilized against dissociation. At 
gy/E/hvABC = 2.1 (Fig. [5b) only the AB + C channel is amplified, and for supercritical 
gy/E/fkVABC — 2.3 (Fig. [5b) atoms are produced in both channels. Control can be 
attained either by varying the coupling strength g or by changing the trap widths. If the 
trap is tightened or loosened so that all trap frequencies are multiplied by a factor w (e.g. 
by variation of the trap field strength) while keeping fixed g, U m and the total number 

(c) (c) 

of particles N, the critical coupling-strength values g\ +BC an d 9ab+c S e ^ multiplied 
roughly by a factor of w 3 ^, because trap level spacing is multiplied by w and the average 
density is multiplied by w 1 ^ (the exact factor in the presence of molecular interactions 
is slightly different, due to variation in the molecular condensate profile). For example, 
for w = 1.2 we have g ( - c \w = 1.2) w 1.2 3 / 4 x g^ c \w = 1) = 1. 1465^^ (w = 1). Retaining 
the dimensionless length and frequency scales of w — 1, we obtain that g^+sc ~ 2.52 
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Figure 6. Channel populations at r = 15 as a function of the widening parameter w. 
The dimensionless time is given in terms of the w = 1 trap period. All parameters are 
the same as in Figs. EJi-lSf- 

(c) 

and g\ B+c ~ 2.29. The actual values from numerical diagonalization turn out to be 
9a+bc = 2.5000 and g^B+c = 2.2680 indicating that the effect of molecular density- 
profile variation is negligible. Similarly, for w = 0.8, the critical values are multiplied 
by 0.8 3 / 4 = 0.8459 to give g^+BC ~ 1-86 and g^B+c ~ 1-69 (in units of the w — 1 trap 
frequency and length), compared to numerical values of 1.8760 and 1.7080, respectively. 
In Figs. [5U-l5f we plot a sequence of dynamical calculations keeping g, U m and N fixed, 
and varying the trap frequency from w — 1.2 (Fig [5H) through w = 1 (Fig[5fe), to 
w = 0.8 (Fig[5f). This sequence corresponds to a realistic experimental scenario where 
only the trap frequency is varied and it clearly shows the dependence of the reaction 
outcome and the channel population ratio, upon the confining potential width. The 
atomic population at r = 15 (in w — 1 dimensionless frequency units) is plotted in Fig. 
[6] as a function w, keeping N, g and U m fixed. It is evident that the ratio between the 
channel populations can be controlled merely by variation of trap width. 

5. Experimental parameters 

Before closing, we assess the relevant experimental parameters required for the 
observation of the predicted supression of stimulated dissociation by tight confinement. 
For the majority of Feshbach-coupled diatomic systems the characteristic atom-molecule 
conversion frequency is much larger than the trap frequency, suggesting they are well- 
within the strong-coupling, nonlinear domain where stimulated gain should be observed. 
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One example is the 23 Na Feshbach resonance at 907 G [TQl tH], where the coupling 
strength is given by g = -JW a [121 W\ with U a = Airh 2 a s /m, k/(2itK) pa 4.6 MHz, 
and a s ~ 60 Bohr. These values give g/2irh ~ 9 x 10 -3 Hz cm 3 / 2 . Given a density of 
n = 10 15 cm -3 , we have g^fnjlTxh pa 284 kHz, far above contemporary trap frequencies 
(we note that comparison of this Feshbach frequency to the atomic interaction frequency 
U a n/27rh ~ 17.5 kHz, justifies the neglect of atom-atom interactions in our calculations, 
even for substantial atomic populations and high densities. For n = 10 13 cm~ 3 , g^fnj^-nh 
drops to 28.4 kHz, whereas U a n/2irh is only 175 Hz). For cj /2vt = 200Hz we obtain 
that N c pa 2 x 10~ 2 , so that stimulated dissociation will be observed for any number of 
molecules in the condensate. It may be that a sufficiently weak Feshbach resonance can 
be found amongst the many resonances of 87 Rb [59] . 

Collective dissociation thresholds, can however be observed in optical stimulated- 
Raman experimental setups [13, [191 EDI EH [22] . In this type of atom- molecule coupling, 
g = Q p Q s /2A 2 p is an effective two-photon Rabi frequency. The one-photon Rabi 
frequencies Q S}P for the pump (inducing a bound-bound transition to an intermediate 
molecular state) and the Stokes (dissociating the intermediate state via a bound-free 
transition) lasers, respectively, are products of laser intensities by overlap integrals 
consisting of electronic transition dipole moments and Franck-Condon factors, and A 2p is 
the two-photon detuning from the intermediate bound state. Thus, the effective coupling 
strength can be controlled with great precision [19] through the adjustment of laser 
parameters. Most appealing are systems of molecules in optical lattices [58 1 IT9 l [211 [22]. 
providing optical coupling to deeply bound internal molecular states, as well as tight 
trap frequencies of u /2ir = 10 — 100 kHz. Currently, these experiments seek to avoid 
collective effects in association by operating in the Mott-insulator regime with unit 
occupation numbers. However, provided that molecular condensates containing roughly 
10 2 — 10 3 particles per site be formed, their dissociation can demonstrate the expected 
amplification threshold. 

For 10 3 87 Rb2 molecules in a 1 kHz trap, the critical value g c = h 2 / {m^J Nl /A 3 ) of 
the interaction strength turns out to be approximately 2%h x 2.2 x 10 -5 Hz cm 3 ' 2 . The 
corresponding average atomic density is of order 10 13 — 10 14 cm -3 and a s ~ 100 bohr for 
87 Rb. Therefore, Un ~ 75 — 750 Hz for a fully dissociated gas, and atomic interactions 
are initially negligible with respect to the trap and coupling frequencies. One possible 
experiment will involve the preparation of a molecular BEC in a trap smaller than the 
resonance healing length, thus arresting dissociation even in the presence of coupling. A 
sudden change either in g or in ujq (which in a lattice can be achieved via switching the 
lattice wavelength or depth) will trigger the stimulated dissociation of the BEC [47j . 

While similar values for the critical parameters apply to the case of heteronuclear 
diatomic molecules, it is more difficult to estimate the pertinent experimental parameters 
for the heteronuclear triatomics described in section HI due to the lack of experimental 
or theoretical values for the required Feshbach resonance strengths. 
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6. Conclusions 

The realization of dilute Bose condensed gases has produced ensembles with 
unprecedented mesoscopic coherence lengths of the order of several \xm. These healing 
lengths are comparable with characteristic confining potential lengths, suggesting a 
wealth of interesting phenomena, such as the stabilization of condensates of attractively 
interacting atoms [52j [53] and the generation of matter-wave solitons. The production 
of molecular BECs [TJ HI [5] opens the way for the utilization of this macroscopic 
coherence to generate a novel type of collective 'superchemistry', in which atoms are 
not treated individually but as one mesoscopic, coherent entity [261 E3 EE]. Previous 
work has shown that bose stimulation effects should dominate atom-molecule dynamics 
in degenerate quantum bose gases. It was found that the resulting exponential gain 
produces extremely heightened selectivity in multichannel processes [28] and can lead 
to pattern formation in molecular dissociation |27j . 

Atom-molecule coupling is characterized by a resonance healing length. If confining 
trap sizes become comparable with this length, the molecular BEC can be stabilized 
against dissociation in the same way an attractively interacting BEC is stabilized 
against collapse [17]. In this work we have used this stabilization effect to show that 
confinement can also be used to control multichannel dissociation. Differences in mass 
and polarizability of the dissociation fragments produce differences in the resonance 
healing lengths and effective trap lengths of the various rearrangement channels, 
resulting in different amplification thresholds, even when the coupling strengths are 
equal. Consequently, the channel population ratio depends critically on the trap width 
and geometry. This effect highlights another novel feature of superchemistry in that 
reaction yields depend on reaction vessels. 
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